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ProG6d6, installation et programme d' ordina-beur pour 
estimer I'effectlf initial d'une population d'acides 
nucleiquefl, notanment par PGR 

5 La presente invention concerne 1 ' estimation de* I'effectif 
initial d'une population d'int6r@t dans un §chantillon 
soumis a succession de reactions d' amplification. 

La pr6sente invention trouve une application avantageuse 
10 mais • aucunement limitative S la dStermination d'une 
quantite initiale d'acides nucleiques dans un 6chantillon 
soumis 4 une reaction de PGR (pour "Polymerase Chain 
Reaction , en temps reel. Une technique de ce type^ dite 
"de quantification par PCR"/ est utilis6e notamment pour 
15 evaluer le nombre de copies d' agents pathogdnes (par 
exemple du Virus de 1' Immunod^f icience Humaine VIH) dans 
un pr^levement de fluides corporels d'un patient, 
typiquement dans le cadre d'un suivi medical, 

20 On se refere a la figure 1 pour decrire bridvement 
1' allure schematique d'une courbe d' amplification par PGR 
en temps reel, avec, en abscisses/ des indices de cycles 
de PGR et/ en ordonneeS/ dans 1 ' exemple represente, des 
quantites de fluorescence ^mises (en unite arbitraire) et 

25 mesur6es ^ chaque cycle de PGR. On indique qu' a chaque 
cycle de PGR, 1' 6chantillon subit des variations de 
temperature- permettant a I'ADN polymerase de proceder ^ 
1' amplification des ac'ides nucleiques et les produits PGR 
correspondants sont . detectes par des molecules 

30 fluorescentes . On obtient, en reportant la quantite de 
fluorescence mesuree Fn en fonction du nombre n de cycle 
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de PGR, une variation du type represente sur la figure 1 
et comprenant au moins : 

- une premiere partie BDF oix les mesures de fluorescence 
se confondent senslblement avec un bruit de fond de 

5 l*appareil de mesure de fluorescence, 

- une seconds partie EXP oix les quantites de fluorescence 
mesurees croissent de fagon senslblement exponentielle, 

- une troisieme partie LIN oix la croissance des quantites 
de fluorescence mesurees est senslblement att6nu6e et 

10 se comporte globalement de fagon senslblement lin6aire, 

et 

- une quatriSme partie PLA oCi les mesures de fluorescence 
atteignent une phase de plateau. 

15 On retiendra en fait que pour les premiers cycles de PGR 
(premidre et seconde parties), la population d'inter§t 
croit senslblement de fa^on exponent ielle, tandis que pour 
les cycles qui suivent (troisieme et quatrieme parties) , 
d'autres phenom6nes entrent en competition avec la 

20 croissance de la population d'interet, cette croissance 
etant alors amortie jusqu'^ la phase de plateau PLA. 

On connait notamment par le document : 

"Mathematics of quantitative kinetic PCR and the 
25 application of standard curves", de R.G.Rutledge et 
C.Gote, dans Nucleic Acids Research , 2003, vol,31, N^'IG, 
une mSthode pour estimer la quantite initiale, inconnue, 
d'acides nucl6iques dans un 6chantillon d'int6ret, par 
PGR. Elle consiste ^ utiliser plusieurs echantillons de 
30 quantites initiales connues en acides nucleiques, dits 
"standards", pour determiner, par interpolation, la 
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10 



15 



20 



30 



quantite initiale d'acides nucleiques pr6sents dans 
I'echantillon d'interfit. 

De fagon g^nerale, plus la quantity initiale en acides 
nucl6iques dans un 6chantillon est 61ev6e, et plus t6t la 
quantity de produit PGR obtenue et, par consequent, la 
fluorescence 6mise deviendront detectables. En se r6ferant 
^ la figure 2 relative ^ I'art anterieur, on comprendra 
que la population initiale dans le standard Stl est plus 
61evee que celle du standard St2, laquelle est plus elev6e 
encore que celle du standard St3, etc, car le cycle Ctl 
pour le standard Stl intervient avant le cycle Ct2 
correspondant pour le standard St2, lequel intervient 
encore avant le cycle Ct3 du standard St3, etc. 



Ainsi, ce cycle Ct, correspondant au cycle auquel la 
mesure de fluorescence atteint un seuil de fluorescence 
THR (comme represents sur la figure 2), pr^fixe 
arbitrairement (typiquement au-dessus du bruit de fond) , 
se presente comme un parametre representatif de I'effectif 
initial No d'une population d' acides nucleiques soumise 
aux cycles de PGR. Gette observation a et6 mise a profit 
dans I'art anterieur cit6 ci-avant pour §tablir une 
d^pendance telle que representee sur la figure 3 entre les 
25 numeros de cycle Gtl, Ct2, Ct3, Ct4, pour plusieurs 
standards d'effectifs initiaux connus, et leur effectif 
initial Nq^, No^, Nq^, Nq*. En reportant ainsi les cycles 
Ctl, Ct2, Gt3, Ct4, etc, en ordonn6es et le logarithme des 
effectifs initiaux NqS No^ Nq^, Nq^ en abscisses, on 
obtient une pente de regression REG. Sur cette pente REG, 
on reporte (fleche Fl en traits pointings) le cycle Ctint 



4 



d6tect6 pour 1' 6chantillon d*interet. Par interpolation 
sur la " pente de regression REG (fleche F2 en traits 
.pointings), on determine ensuite I'effectif initial Nq^"^ 
de la population dans 1' echantillon d'interSt. 

5 

Cette mSthode, mSme si elle s'est largement repandue, 
pr6sente neanmoins quelques inconvenients . 

Tout d'abord, elle necessite 1 ' utilisation de plusieurs 
10 6chantillons standard d^effectifs initiaux respectifs 
connus . 

Ensuite, la methode depend du juge de 1' utilisateur car la 
valeur du seuil de fluorescence qui est choisie par 

15 1' utilisateur influence directement la valeur des cycles 
Ct des courbes d' amplification et, par consequent, la 
valeur estim^e de I'effectif initial dans 1' Echantillon 
d'interet. La valeur du seuil a egalement un impact sur la 
precision du resultat car la precision est g^neralement 

20 meilleure si le seuil est choisi dans la phase EXP de 
croissance exponentielle des courbes d' amplification . 
Toutefois, en pratique, il est difficile pour 
1' utilisateur de savoir si le niveau seuil de fluorescence 
fixe THR correspond bien a la phase exponentielle des 

25 courbes, et ce, pour tous les Echantillons (les 
6chantillons standard et 1 ' Echantillon d'interet). 

Enfin, la mEthode sous-entend, sans aucune verification, 
que la population a le meme rendement d' amplification dans 
30 1' Echantillon d'interet et dans tous les echantillons 
standard. Ainsi, si typiquement 1' echantillon d'intEret 



Copy provided by USPTO from the IPW Image Database on 01/19/2006 



contient des inhibiteurs de PGR, son resultat sera 
- faussement abalssS. 

On retiendra done que la technique de I'art ant§rieur 
5 depend du seuil de fluorescence THR dSfini par 
1 'utilisateur . La valeur choisie influence les valeurs des 
cycles Ct * et par consequent la determination de • la 
quantity initiale dans 1 ' 6chantillon d'int6r§t, C'est 
I'une des raisons pour lesquelles de nombreux travaux ont 
10 cherche receiranent a automatiser et ^ fiabiliser la 
detection du cycle Ct. 

La presents invention vient ameliorer la situation en 
proposant une approche compietement differente. 

15 

Elle vise tout d'abord un proc6d§, ce proc6d§ etant mis en 
ceuvre par des moyens inf ormatiques pour quantifier, de 
fa^on absolue et/ou relative, une population initiale 
d'acides nucleiques dans un echantillon d'interet. Get 

20 echantillon est soumis a une succession d ' applications 
d'une reaction d' amplification de la population d'interet. 
De fagon tres generale, cette amplification peut etre 
men6e par cycles successifs de PGR, mais aussi par toute 
autre technique d' amplification. On retiendra surtout que 

25 cette amplification doit simplement etre definie par un 
rendement de reaction, comme on le verra plus loin. Au 
cours de ces amplifications successives, on releve des 
mesures experimentales representatives d'un effectif 
courant de la population au moins dans 1 ' echantillon 

30 d'int6r§t. On comprendra que I'on peut relever une ou 
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plusieurs mesures apres ou pendant chaque reaction 
d' amplification^ sans perte de g6n§ralite. 

Selon une definition actuellement pr§f er6e de 1 • invention^ 
5 le precede au sens de 1' invention comporte les 6tapes 
suivantes : 

a) pr6voir un modfele de rendement de la reaction 
d' amplification en fonction de la succession 
d' amplifications, ce modele comportant : 

10 o une phase . sensiblement constante pour une premiere 

partiedes applications de la reaction, et 
o une phase non constante pour une seconde partie 
des applications de la reaction, 
les premiere et seconde parties 6tant s6par6es par une 
15 region de basculement du rendement entre les phases 

constante et non constante, incluant un indice 
d' application correspondant sensiblement au 

basculement^ 

b) exprimer, en utilisant le modele de rendement, une 
20 relation faisant intervenir au moins 1 • indice de 

basculement et un parametre representant I'effectif 
initial dans I'echantillon d'interet, 

c) par comparaison avec les mesures experimentales r 
determiner au moins 1' indice de basculement, et, dans 

25 une etape d) ult^rieure ou conjointe, en deduire 

I'effectif initial dans 1 ' 6chantillon d*int§r§t. * 

D'autres avantages et caracteristiques de 1 " invention 
apparaltront A la lecture de la description d6taill6e du 
30 mode de realisation donnS ci-apres d titre d'exemple, en 
s'appuyant sur les graphiques annexes et sur lesquels : 
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- outre la figure 1, relative S I'art ant^rieur, 
repr6sentant une variation de la quantity de 
fluorescence mesurSe en fonction d'un noxnbre de cycles 
de PGR, decrite ci-avant ; 

5 - outre la figure 2, relative k I'art ant^rieur, 
representant les quantites de fluorescence §mises et 
- croissantes en fonction du nombre de cycles de- PCR^* et 
d6crite ci-avant ; 

- outre la figure 3 representant • une mSthode 
10 d* interpolation pour dfeterminer la quantity initiale de 

la population d'int6ret dans 1 ' 6chantillon d'int6r§t, 
selon une m6thode connue de I'art anterieur, et d6crite 
ci-avant ; 

- la figure 4A represente une variation schematique des 
15 mesures experimentales pr6citees en fonction de la 

succession d ' amplifications que I'on applique k 
1 ' 6chantillon d'interet, 

- la figure 4B represente une variation sch6matique du 
rendement de la reaction d' amplification, obtenue a 

20 partir des mesures experimentales, en fonction de la 

succession d ' amplifications appliqu6es ^ 1 ' 6chantillon 
d ' interet , 

- la figure 5 represente une regression entre des indices 
de basculement du rendement entre la phase constante et 

25 la phase non constante et le logarithme des effectifs 

initiaux pour des echantillons standard et pour 
1' echantillon d'interet, pour la mise en oeuvre d'un 
premier mode de realisation, 

- la figure 6A represente une variation typique de la 
30 quantity de fluorescence mesuree et ajustfie en tenant 
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compte d'un bruit de fond propre aux mesures, en 
fonction du nombre de cycles n de la reaction PGR, 

- la figure 6B repr6sente une variation de I'efficacite 
de la PGR de la figure 6A, en fonction du nombre de 

5 cycles n, 

- la figure 7 represente une comparaison de la variation 
exp&rimentale de la fluorescence 6mise selon la figure 
6A et le rfesultat obtenu par 1 ' application du meddle de 
la fluorescence emise obtenu par integration d*un 

10 module d' ef f icacite, selon un second mode de 

realisation, 

- la figure 8 repr§sente une comparaison de la variation 
de I'efficacitS de la figure 6B avec 1 ' application d'un 
module d'ef f icacite, d^duit de la mod61isation de la 

15 fluorescence emise de la figure 1, 

- la figure 9 repr§sente 1 ' organigrarame schematisant les 
principales etape's du proc6d6 selon le mode de 
realisation propose pour la presente invention, et 

- la figure 10 represente une installation pour 
20 quantifier la population initiale d'un echantillon 

d'interet • 

On s'appuie sur les figures 4A et 4B pour decrire 
brievement quelques principes de 1' invention illustrant 
25 les caracteristiques du precede ci-avant. 

Tout d'abord on indique que la figure 4 A represente une 
succession de mesures 'exp^rimentales Fn representatives de 
I'effectif courant d'une population d'interet qui subit 
30 progressivement une succession de reactions 

d" amplification, chaque reaction etant indexee par un 
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indice n. Cette succession de reactions correspondr dans 
I'exemple decrit ici ^ titre non limitatif^ ^ une 
succession de cycles de PGR, Les mesures expferimentales Fn 
correspondent dans cet exemple ^ titre non limitatif a des 
5 mesures de quantit§s de fluorescence a chague cycle de 
PGR- Ainsi, dans le cadre d'une m6thode de quantification 
qui combine la reaction PGR et la fluorescence emise par 
1 ' fechantillon d'inter^t, on introduit des reactifs 
fluorescents dans 1 ' 6chantillon, de sorte que la 

10 fluorescence 6inise, au cours d'un cycle de PGR, est 
proportionnelle d I'effectif d'acides riucl6iques dans cet 
echantillon. N6anmoins, cette realisation admet de 
nombreuses variantes, Deja, il peut etre pref6r6 de mener 
plusieurs mesures ou aucune mesure pour certains cycles de 

15 PGR. Par ailleurs, plus g^neralement , la methode de mesure 
peut utiliser d'autres techniques que la fluorescence, 
meme si cette derni^re est souvent celle employee dans la 
quantification par PGR. Enfin, on indique que d'autres 
techniques d' amplification peuvent etre prevues dans le 

20 cadre de la presente invention, dds lors qu'il est 
possible de suivre la variation du rendement de la 
reaction correspondant ^ cette amplification. Neanmoins, 
comme 1* exemple decrit ci-apres traite pref §rentiellement 
de cycles de PGR, on parlera efflcacite de la PCR'\ 

25 not6e En pour chaque cycle de PGR d' indice n, pour 
designer le rendement de la reaction d ' amplification. 

Comme indiqu6 ci-avant en reference ^ la figure 1, on 
retrouve sur la figure 4A principalement deux regions ot : 
30 - pendant les premiers cycles de PGR (partie EXP) , la 
population croit sensiblement de fagon exponentielle. 
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- tandis que pour les cycles qui suivent (partie LIN et 
PLA) , d'autres phenomenes entrent en competition avec 
la croissance de la population d'int6r§t, cette 
croissance 6tant alors amortie. 

On forinule alors les hypotheses suivantes : 

- le rendement de la reaction En est relativement constant 
durant les premiers cycles sur la partie EXP, et 

- a partir d'un certain nombre de cycles effectu6s, le 
rendement En de la r6action decrolt sur la partie LIN et 
PLA • 

Cette decroissance du rendement peut avoir diverses 
explications, notamment une degradation et/ou une carence 
des reactifs de PGR (DNA polymerase, dNTPs, amorces, ou 
autres) et/ou une inhibition par les produits form6s eux- 
memes - 

On fait ainsi I'hypothese ici que le rendement est d'abord 
constant et decroit ensuite, Toutefois, on comprendra que 
1' invention s* applique de fagon plus . gen^rale dans le 
cadre d'un rendement : 

- d'abor.d constant, ce qui correspond A une situation 
normale de croissance par amplification, 

- puis non constant (decroissant ou croissant) , ce qui 
correspond a une situation sensibleirient anormale. 

On a constate que dans le cadre de reactions 
d' amplification de quantite d'acides nucleiques, le 
rendement bascule souvent d'une phase constante a une 
30 phase non constante. Cette observation a ete mise S profit 
au sens de 1' invention pour en d6duire la quantite 



20 



25 
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initials d'acides nucl6iques, coinme on le verra en detail 
ci-apr6s. Auparavant, on indique simplement que le 
rendement peut encore basculer d'une phase non constante 
pour les premiers cycles a une phase constante ensuite. La 
5 presente invention s' applique tout aussi bien a ce cas, De 
maniSre g6n6rale, on retiendra done que I'on detecte, au 
sens de 1' invention, un basculeraent du rendement; d'une 
phase constante k une phase non constante. 

10 On rappelle que le but recherche est de donner I'effectif 
initial de la population soumise S amplification. D6ja, en 
se r6f6rant a la figure 4A, on comprendra que la mesure Fo 

, . representative de cet effectif initial, pratiquement 
confondue avec le bruit de fond des mesures BDF ne permet 

15 pas directement . d'acceder, k elle seule, a 1' effectif 
initial. Dans I'art anterieur, on a recherche a quantifier 
cet effectif initial en exploitant la phase exponentielle, 
done typiquement a la sortie du bruit de fond. On y 
determine alors le cycle seuil Ct (correspondant au point 

20 AA pour "art anterleur") de la figure 4A. Or, comme on I'a 
vu ci-avant, les mesures sont souvent bruit6es dans cette 
region et il est difficile de determiner avec precision un 
cycle seuil Ct de sortie du bruit de fond. 

25 Dans une toute autre approche, la presente invention 
exploite plutdt la quasi-totalit6 des points de la courbe 
d' amplification pour determiner pr6cis6ment une region BAS 
oCi le rendement bascule entre une phase constante et une 
phas^ non constante, typiquement ici entre la phase 

30 • exponentielle EXP et la phase lin^aire LIN. Deja, on 
constatera que les mesures sont logiqueraent moins bruitees 
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dans cette region BAS que dans la region de sortie du 
bruit de fond puisque la region BAS intervient pour des 
cycles ultSrieurs. De plus, notamment en raison des 
propri6t6s mathSmatiques li6es au rendement, on verra plus 
5 loin que le nombre de standards qui peuvent Stre utilises 
pour la quantification de I'effectif initial dans la 
•population - d-'interet est, de ' fagon* tres -avantageuse; 
r6duit par rapport au nombre de standards utilises dans la 
quantification de l*art anterieur. • 

10 

On d6crit maintenant briSvement la relation qui peut lier 
la region de basculement BAS a I'effectif initial de la 
population d'int^ret. On rappelle au prealable que le 
rendement d'une reaction d' amplification est donn6 par : 

15 Nn+l = Nn + En.Nn, OU : 

- Nn est I'effectif de la population d'interet apres une 
amplification d'indice n dans une succession 
d ' amplifications, 

- Nn+i est I'effectif de la population d'interet apres une. 
20 amplification suivante, d'indice n+1, dans la 

succession d' amplifications precitee, et 

- En est le rendement de la reaction d ' amplification 
d'indice n dans la succession d • amplifications 
precit6e. 

25 En reformulant cette relation, par recurrence, on 
obtient : 

N„+i = {1 + En) (1 + En-l) (1 + En-2) d + Eq) Nq, 

OU No est I'effectif initial de la population d'int^r^t. 
Tant que le rendement En est constant, on comprendra que 
30 cette derniSre relation peut s'ecrire simplement : 



f 
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Nn+i - (1 + Eo)"-*-^ . No (Oil I'indice n+1 n'a pas 
encore atteint la region de basculement BAS) • 
En effet, le rendement 6tant alors constant pour ces 
premiers cycles, on a : 

5 En « En-i = En-2 ^ ... = Eq, 

oCi Eo est la valeur de la phase constante du rendement. 
Neanmoins/ lorsque I'indice n+1 depasse la r6gion de 
basculement BAS, la relation devient : 

Nn+i = No X (1 + Eo)^EEP X fonation (Cebp^ oCi : 

10 - (Cbep " 1) est le dernier indice de la reaction 
d' amplification au cours duquel le rendement est encore 
constant (on comprendra done que I'indice Ceep 
represente lui-meme I'indice de basculement, proprement 
dit, entre la phase de exponentielle k la phase 
15 lineaire) , 

- le terme f onct Ion (Ceep r n+1) est une fonction 
particulidre qui caract^rise la phase non constante du 
rendement et qui depend au moins de I'indice de 
basculement Ceep et de I'indice courant de 
20 1 ' amplification n+1. 

On voit ainsi comment peuvent §tre lies I'indice de 
basculement Ceep et I'effectif initial No de la population 
d'interet. On comprendra qu'a ce stade, on a d6ja realise 
25 les etapes a) et b) du precede tel que d^fini ci-avant. 

D6ja, un premier mode de realisation consiste ^ determiner 
experimentalement I'indice de basculement Ceep' et ^ le 
corr61er a I'effectif initial, par regression, en 
30 utilisant plusieurs echantillons standard qui subissent le 
meme traiteraent d ' amplification que 1 ' 6chantillon 
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d'int6r§t,' Dans ce cas, on comprendra que les Stapes b) et 
c) du proc6de tel que d6fini ci^avant sont simplement 
inversfees puisque I'on determine exp6rimentalement, 
d'abord, I'indice de basculement Ceep (etape c) ) et que 
5 I'on determine ensuite la relation entre I'indice Cbep et 
I'effectif initial No (etape b) ) pour obtenir finalement 
I'effectif -initial No (etape d) ) . ... 

Avant de decrire en detail toutes ces Stapes au sens de ce 
10 premier mode de realisation, on dScrit un precede de 
determination de I'indice Cebp/ a partir des mesures 
expSrimentales . En particulier, on comprendra que cette 
manidre de determiner experimentalement I'indice Ceep peut 
Stre appliquSe A un autre mode de realisation, different 
15 du premier mode de realisation prfecite, 

En reprenant la relation entre I'efficacitS En A un cycle 
donne n et I'effectif courant dans 1 ' echantillon d'intSret 
au mSme cycle Nn et a un cycle ultSrieur Nn+i, on peut 
20 exprimer I'efficacite de 1 ' amplification comme suit : 

En = Nn+l/Nn " 1 

Dans certains cas notamment ou il n'est pas necessaire de 
tenir compte du bruit de fond BDF dans les mesures, on 

25 peut en premiere approximation considSrer que les mesures 
sont sensiblement proportionnelles a I'effectif courant 
dans 1 • echantillon d'interet. Toutefois, en pratique, on 
tiendra le plus souvent compte des derives des mesures et, 
dans une realisation qui sera decrite maintenant, on 

30 determine des . mesures experimentales corrigSes, notees 
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F'nr a partir des mesures directes Fn telles que 
reprfesentSes sur la figure 4A, 

Pr§f§rentiellement, on applique une §tape pr6alable de 
5 traitement des mesures exp6rimentales Fn consistant S 
soustraire le bruit de fond BDF et a introduire ensuite 
une compensation pour tenir compte d'une mesure non nulle 
e representative de I'effectif initial. Dans I'exemple 
repr6sent6 sur la figure 4A, la variation du bruit de fond 

10 BDF en fonction de I'indice n peut Stre representee par 
une fonction lin6aire car les essais ont raontr6 qu'un 
modele lineaire est satisfaisant pour des mesures de 
fluorescence en PGR. Neanmoins, un module de variation 
exponentielle pourra etre pr6f§re dans certains cas , On 

IS applique en tout etat de cause le module qui respecte le 
mieux la variation du bruit de fond BDF donnee typiquement 
par les premiers points de mesure. Ensuite, on soustrait 
le modele BDF de variation de bruit de fond choisi a 
toutes les valeurs de mesures exp6rimentales Fn. En 

20 appliquant cette etape, on comprendra que la mesure de 
fluorescence theorique Fo est ramenee a une valeur de 
mesure nulle, ce qui correspond ^ un effectif initial No 
nul, ce qui n*a aucune realite physique. Par consequent, 
on applique avantageusement une compensation de cette 

25 correction comme suit : 

F'n = Fn - BDF + ou : 

- le terme F'n correspond a une mesure corrigee pour un 
indice courant n, 

- le terme Fn correspond a la mesure experimentale brute a 
30 cet indice courant n. 
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- le terme BDF correspond A la valeur du bruit de fond 
mod61is6 pour I'indice courant n^ et 

- e est le terme de compensation correspondant , supposS 
constant dans 1' example d6crit et representant 

5 directement en fait I'effectif initial Nq. 

Ces 6tapes de correction par rapport au bruit de fond, 
meme si elles sont trds avantageuses pour la determination 
de I'indice de basculement Cebp/ peuvent §tre appliquees 
par ailleurs S toute determination et quantification de 
I'effectif initial No d6s lors qu'un bruit de fond est 
susceptible de fausser la mesure de cet effectif Nq. A ce 
titre, ces etapes peuvent faire I'.objet d'une protection 
s^paree, le cas 6ch6ant, 

Les mesures corrigees F'n ainsi obtenues sont 
avantageusement proportionnelles aux effectifs courants Nn 
dans 1 'echantillon d'intferet, de sorte que le rendement En 
peut maintenant s ' exprimer directement en fonction des 
valeurs de mesure (corrigees comme indiqu6 ci-avant) , par 
la relation : 

En - F'n+l/F'n " 1 

Ainsi, ^ partir des mesures experimentales Fn de la figure 
4A, on obtient les mesures corrigees F'n a partir 
desquelles on determine ensuite la variation de 
I'efficacite En en fonction de la succession d' indices n, 
telle que representee sur la figure 4B. 

Eh bref , on exprime les mesures experimentales sous la 
forme d'une variation experimentale du rendement En telle 
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que representee sur la figure 4B, en fonction de la 
succession d' applications n. On obtient alors une 
variation experimentale du rendement, comportant : 

- une premiere region sensiblement bruit§e pour des 
5 indices d' applications n faibles (sensiblement avant le 

cycle CG dans I'exemple de la figure 4B) , 

- suivie d'une seconde region moins bruitfee pour des 
indices d' applications plus eleves (au . moins aprSs la 
region de basculement BAS) . 

10 

Au moins dans le cas le plus courant en amplification par 
PGR et en mesure par fluorescence, la phase non constante 
du rendement est decroissante et correspond a cette 
seconde region peu bruitee (comme repr6sent6 sur la figure 
15 4B) . Notanunent pour eliminer des points de mesure qui 
risqueraient de faiisser les r6sultats lorsque I'on choisit 
un module ^ appliquer k la variation du rendement : 

- on estime une valeur grossi^re de la phase constante du 
rendement Eq^ et 

20 - notamment pour la recherche de 1 ' indice de basculement 
CsEPr on ignore une partie au moins des mesures dans la 
seconde region moins bruit§e et dont le rendement 
estime est inferieur a une valeur seuil, telle que, par 
exemple, une fraction de la phase constante Eq. 

25 

Ces points NEG (figure 4B) qui sont elimines sont 
typiquement ceux correspondant a des indices 
d' amplification n tres eleves et qui pourraient ne plus 
respecter le modele qui a ete choisi pour I'efficacite 
30 sensiblement autour de la region de basculement BAS. Pour 
les 61iminer, on ^valuera par exemple une moyenne de la 
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phase constante du rendement Eo typiquement pour les 
premiers indices n. Ensuite^ on choisira une valeur seuil 
correspondant ^ une fraction," par exemple 10%, de la 
moyenne trouv6e pour la phase constante Eq. Ensuite, en 

5 partant des indices n les plus eleves, on elimine tous les 
points de mesure NEG dont la mesure de rendement est 
inf6rieure ou 6gale ^ cette valeur* seuil.* Cette 6tape, 
tr6s avantageuse pour la detection de I'indice CsBPr peut 
n6anmoins etre appliquee a toute determination basee sur 

10 le rendement En et, le cas echeant, peut faire I'objet 
d'une protection s6par6e. 

Dans le cas d'une phase non constante du rendement qui est 
dScroissante et qui suit la phase constante, coirane 

15 repr6sent6 sur la figure 4B, on rep§re la region de 
basculement BAS, dans le sens des indices n d6croissants a 
partir de la seconde region moins bruit §e, en dStectant un 
indice grossier CG dont le rendement depasse une valeur 
predetermines. Ainsi, en se referant a la figure 4B, on va 

20 dans le sens des indices n decroissants en remontant vers 
la region de basculement BAS et .on evalue le rendement 
associe a cet indice. Pour le premier point de mesure dont 
le rendement depasse signif icativement la valeur 
predeterminee precit^e, on considere que I'indice grossier 

25 precite CG est detects et correspond a I'indice de ce 
point de mesure . 

Comme on le verra dans un mode de realisation ulterieur, 
on peut modeiiser, pour chaque point de mesure, la 
30 variation de son rendement comme si ce point correspondait 
ef fectivement a I'indipe de basculement Ceep- Dans cette 
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realisation, on estime alors sa phase constante Eo de 
rendement et si cette valeur estim6e d^passe la valeur 
pred6termin6e pr6cit6e, on considdre ce point coinme 
correspondant d I'indice grossier CG, 

5 

De fagon generale, un rendement maximum a une valeur de 1 
de soxte- que I'on peut choisix la. valeur pr§d6termin6e 
pr6citee comme etant egale A 1, Neanmoins, ce choix admet 
des varlantes et par exemple, il peut §tre prevu de fixer 
10 cette valeur pr6d6termin§e comme correspondant k la 
moyenne du rendement Eq 6valu6e sur les premiers cycles de 
reaction . 

Ensuite, on affine 1' estimation de la valeur, qui peut 
15 avantageusement etre f ractionnaire, de I'indice 
d' application Ceep dans la region de basculement, dans le 
sens des indices d' applications croissants, a partir de 
I'indice grossier CG, en detectant un indice d' application 
d* amplification dont le rendement vaut environ la valeur 
pr6d6terminee pr6citee, Ainsi, en se r6ferant a nouveau a 
la figure 4B, pour af finer la recherche de I'indice de 
basculement Ceep apres la determination de I'indice 
grossier CG, on redescend a partir de I'indice grossier CG 
dans le sens des indices n croissants, par pas plus petit 
qu'un indice entier et on determine, par exemple par 
interpolation, I'abscisse pour laquelle on franchit la 
valeur pr6d6termin6e . Typiquement, tant que la valeur de 
phase constante reste superieure a 1, on poursuit la 
recherche dans le sens des n croissants et I'on determine 
I'indice (Ceep - 1) qui precede le basculement des que la 
valeur constante Eq est egale ou est tres proche de 1. 



20 



25 



30 



Copy provided by USPTO from the IFW Image Database on 01/19/2006 



10 



15 



20 



2S 



30 



20 

de 



C'est la raison pour laquelle on choisira un pas _ 
recherche correspondant d une fraction d'indice, par 
exeraple 10% d'un cycle n. 

Dans le premier mode de realisation pr6cit6, on pr6voit 
une pluralite d • 6chantillons standard d'effectifs initiaux 
respectifs nconnus, auxquets on applique la successron 
d' amplifications, sensiblement dans les memes conditions 
que pour 1 • ^chantillon d'int§r§t. On determine aussi leurs 
indices respectifs de basculement selon les 6tapes a), b) 
et c) decrites precedemment . A l'6tape d) : 
- on etablit une dependance entre les effectifs initiaux 
des 6chantillons standard No='= et leurs indices Cebp", 
et, apres la determination de I'indice pour 
l'6chantillon d'int6r§t Ceep, on determine I'effectif 
initial de la population d'int6rSt No par interpolation 
sur cette dependance. 

Ainsi, en se rSf^rant a la figure 5, on peut 6tablir une 
dependance entre le cycle de basculement W,Ceep2 des 
standards et leur concentration initiale No\No2 (en fait 
leur logarithme), par exemple par regression. En mesurant 
I'indice de basculement Ceep^"'^ pour 1 ' echantillon d'interSt 
et en reportant sa valeur sur la pente de regression de la 
figure 5, on obtient par interpolation la concentration 
initiale No^"'= dans 1 ' echantillon d'int6r§t. 

Ce premier mode de realisation s' apparente done 
sensiblement k celui de l-art anterieur decrit en 
reference A la figure 3. Cependant, on gardera k 1' esprit 
que I'indice de basculement Cbbp sur lequel se base ce 
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premier mode de realisation ne correspond aucunement au 
cycle Ct de I'art ant6rieur. 

Dans une approche sensiblement differente de ce premier 
5 mode de realisation : 

- a I'^tape b) , on exprime en fait, en utilisant le 
module *de rendement, une" variation parametrfee en 
fonction de la succession d ' amplifications, cette 
variation mettant en jeu au moins un parametre 

10 repr6sentant I'indice de basculement CEEPr et 

- & I'etape c) , par comparaison avec les mesures 
exp6rimentales, on determine au moins ce parametre 
repr6sentant I'indice de basculement Cbep. 

15 Dans un second mode de realisation, cette variation 
parametree est representative de I'effectif courant Nn 
dans 1 'echantillon d'interet, 

Typiquement, cette variation parametree peut etre tiree de 

20 1' expression du type donne ci~avant : 

Nn+i = No X (1 + Eo)^EEP X fonction (CEEPf n-f-l) 
Ainsi, cette variation met en jeu, outre un parametre 
repr^sentant I'indice de basculement Ceep/ un parametre 
representatif de I'effectif initial dans 1 • echantillon 

25 d' inter et No- 

Ensuite, aux etapes c) et d) de ce second mode de 
realisation, on determine sensiblement con jointement ces 
deux parametres Cbep et No- 

30 
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Auparavant, S l'6tape a), il faut determiner un module 
pour la fonction precitee fonction (CsBPf n+1) . 

Le plus souvent, pour une quantification par PGR, on 
5 choisira un modele pour la phase non constante du 
rendement correspondant k une exponentielle d6croissante 
comportant un param^tre de" decroissance -(J -qui sera- ptus 
d6taill6 ci-apr6s. Ce paramdtre de decroissance p est 
alors d6termin6, k l"6tape c) , au moins avec I'indice de 
10 basculement Cebp, par comparaison avec les mesures 
experimentales • 

Ainsi^ dans ce second mode de realisation, une fois le 
modele du rendement En choisi, on 1' applique k 
15 1' expression generale de I'effectif courant Nn donnee par 
la relation ci-avant. On obtient alors un module de la 
variation de I'effectif courant Nn- 

Toutefois, ^ moins que les mesures experimentales ne 
20 donnent directement la valeur de I'effectif courant Nn (ce 
qui est rarement le cas en pratique^ actuellement ) , il 
convient de modeliser ensuite les mesures experimentales 
Fn/ elles-m§mes, en gardant a 1' esprit la soustraction du 
bruit de fond et la compensation consecutive e decrites 
25 ci-avant. 

Ainsi, dans une realisation actuellement preferee, la 
variation parametree precitee : 

- est representative des mesures experimentales, et 
30 - comporte un parametre correspondant a une valeur de 
mesure Fo representative de I'effectif initial. 



Copy provided by USPTO from the IFW Image Database on 01/19/2006 



23 

On determine ensuite la valeur de mesure de I'effectif 
initial Fq par comparaison de cette variation parametr6e 
Fn avec les mesures experimentales . 

5 En fait, pour mener cette comparaison, on pourra par 
exemple ajuster les parametres Fo, Eo, Ceep et le parametre 
de d^croissance p-dans le modSle des mesures Fn -par 
correlations statistiques (typiquement la methode des 
moindres carr§s) qui s'appuient sur les mesures 

10 experimentales brutes. Un exemple de realisation est 
d6taill6 ci-apr6s , 

On obtient d'abord une variation de la quantity de 
fluorescence mesuree et ajustee, en fonctlon du nombre de 

15 cycles de PGR appliques, telle que representee par exemple 
sur la figure 6A. Cette figure montre la courbe 
d' amplification d'un echantillon d'inter§t contenant des 
acides nucleiques, ici un fragment d'ADN de quantity 
initiale de 100000 copies, marqu6 par 1 ' intercalant 

20 SYBRGREEN lors de la reaction de PGR que I'on effectue sur 
I'app'areil I-Cycler IQ® de la society BIO-RAD®. 

Dans 1' exemple decrit, on comprendra que la reaction 
d' amplification est une reaction de PGR en temps reel, Les 
25 mesures experimentales representent des quantites de 
fluorescence emises . 



On note ci-aprds Fn 
pr6alablement ajustee pour 
30 ci-avant. On note aussi 
theorique avant le premier 



la fluorescence au cycle n 
le bruit de fond, comme indiqu6 
Fp la fluorescence initiale 
cycle. On note En I'efficacite 
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de la PGR au cycle n. On designe par N le nombre total de 
cycles effectu6s durant la reaction de PGR. 

Par hypoth^se, la^ fluorescence mesuree ^ chaque cycle n de 
la reaction de PGR est definie par : 

FnH « Fn(l-^EJ pourtoutn6{0.1,2,...,N-l} (1) 
avec 0 ^ En ^ 1 . 

On calcule I'efficacite de la reaction a chaque cycle n, 
conune suit : 

E„ = - 1 « ^ - 1 pour tout n€{l,2,...,N.l} (2) 



On remarquera que l"6quation (1) est suppos6e vraie pour 
n=0. Toutefois, par definition, la fluorescence initiale 
Fo est inconnue- L'efficacit6 au premier cycle Eo ne peut 
done pas Stre calcul6e directement par la formule (2). 

On a repr6sent6 sur la figure 6B I'efficacite de la 
reaction PGR approxim§e par la formule (2) et A partir de 
la variation de la fluorescence ajustee de la figure 
en fonction du nombre de cycles n. 

On formule pref erentiellement les hypotheses suivantes : 

- I'efficacite de la reaction est relativement constante 
durant les premiers cycles, et 

- a partir d'un certain nombre de cycles effectues, 
I'efficacite de la reaction decroit. 

La figure 6B vient conforter la seconde hypothese, puisque 
I'on peut constater la decroissance de I'efficacite a 



25 

partir du cycle n=17 . En revanche, les efficacit6s des 
cycles 1 a 16 sont fortement bruitees ce qui rend 
difficile la verification graphique • de la premiere 
hypothese • 
5 • 

Toutefois, on affecte pr6f erentiellement a la variation de 
1 ' ef-f icaci"t6 "un moddle du * type* incluant : 

- une premidre phase, constante, entre le premier cycle 
de PGR et le cycle (Cbep-1) qui precede le cycle de 

10 basculement not6 Ceep* et 

- une seconde phase, dScroissante, pour les cycles 
sup6rieurs ou 6gaux au cycle (Ceep-1) - 

Le cycle (Cebp-1) represente done le dernier cycle, qui' 
15 peut dtre f ractionnaire, auquel I'efficacite est encore 
constante . 

On propose alors de mod61iser l'efficacit6 de la reaction 
de la fagon suivante : 

Eo pour O^n ^{C^-l) 

" Ml + Eor^-«"^-^^^> - 1 pour (Chp - 1) ^ n ^ (N-l) ^ ^ ^ 

oil Eo et p sont des paramfetres r6els qui sont estim6s 
grace S la courbe d' amplification de la figure 6A, ou 
encore grSce a la courbe d'efficacit6 de la figure 6B, 
comme on le verra ci-aprSs, 



25 



En variante, on indique qu'un autre choix, par exemple 
parmi les modeles ci-apres Fl i F3, peut etre pr6f6r6, 
notamment en fonction du type d'acide nucl6ique que I'on 
souhaite quantifier : 
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Fl : E„ -exp(.p(n.CEEp+l)).l 
F2 : E„ - exp (-hCii-Ceep+I)") - 1 
F3 : E„ = Of - exp i'liin^C^+\y) - 1 



On estime pr6f6rentiellement plusieurs jeux de paramStres 
^ I'etape c) pour plusieurs candidats de cycle de 
5' basculemenr CEEPf et I'on choi'sit le cycle candidat iillnimal 
dont les paramdtres associSs maximisent les correlations 
statistiques qui peuvent §tre menses a I'etape c) , er> tant. 
que cycle de basculement Ceep. 

10 Comme on I'a vu pr6c6demment, 1' expression (1) ci-avant 
peut encore s ' ecrire sous la forme : 

Fn = Fo flO+Ek) Pourn e {1,V„N} (4) . 

Ainsi, en introduisant I'expression de l'efficacit6 (3) 
15 dans la formule (4), on obtient un nouveau modele a quatre 
parametres {Fo, Eq, P, Ceep) pour la fluorescence 6mise 
ajust6e Fn : 



Fo(l+Eo)" pour l:gn £ C^^ 

FqO+Eo) ^ «pO)-« pour Ceep ^ n < N 



20 L'effectif initial dans la population d'interet No/ 
I'efficacite Eo de la reaction ^ n=0, le param6tre p et le 
cycle de basculement Ceep/ sont lvalues repetitivement pour 
plusieurs valeurs de cycles dans la region de basculement 
BAS, pour trouver un maximum de correlations statistiques 

25 atteint pour une valeur de cycle minimale egale au cycle 
de basculement Cbep- 
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Dans ce second mode de realisation, on pr6f6re modfeliser 
la variation des quantitis de fluorescence mesur6es et 
ajustSes, en fonction des cycles, en se basant sur la 
5 modelisation de la variation de 1 ' ef f icacite, et mener 
ensuite les correlations directement sur les quantites de 
fluorescence mesur^es et ajustees". 

On remarquera qu'en ajustant les mesures de fluorescence 
10 Smises pour le bruit de fond, on cr6e aussi, de fagon 
factice, un ajustement de la fluorescence initiale Fq- 
Ainsi, estimer les parametres du modele d*efficacit6 a 
partir de mesures d'efficacit6 qui sont deduites des 
mesures de fluorescence ajustees, est une source d'erreurs 
15 suppl6mentaires et il peut etre pref §r§ de proceder . en 
deux etapes, comme on le verra dans le troisi^me mode de 
realisation qui sera decrit plus loin. 

Toutefois, le second mode de realisation propos6 est plus 
20 simple et s'adapte bien aux quantifications par PGR 
utilisant des mesures de fluorescence. II se base sur les 
mesures. reelles de fluorescence F' n/ qui correspondent aux 
mesures de fluorescence ajustees pour la derive du bruit 
de fond avec une compensation e menee sur ces mesures. Le 
25 bruit de fond ayant et6 soustrait, on retrouve une 
relation du type : 

F'n = Fn + 8. 

ou e est une quantite qui peut dependre ou non du cycle n. 
Elle est choisie ici constante, pref erentiellement . 

30 
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Dans ce cas, I'efficacite mesuree et "ajustee" aussi, 
notee E'n au cycle, n est d^finie par : 

E'n - ^ - 1 - - 1 pour tout n e {l,2,....N.l} ( 7 ) 



5 Le modele de la relation (5) ci-avant devient : 
_ J F'oO'*"E\,)" -8 pour l£n ^ CgEp 

^ F', - «p(pv 1 . e pour C^^, Sn^N 

Dans ce cas, les valeurs d'efficacite E'n sont approxim^es 
exp6rimentalement ^ partir des mesures, de fagon a pouvoir 
10 fixer un seuil d'efficacite minimal tolere dans la phase 
d'efficacit6 decroissante . On determine ainsi un cycle 
seuil au-del^ duquel les mesures de fluorescence ajustees 
(points NEG figure 4B) ne seront pas exploitees pour mener 
la modelisation. Typiquement, le cycle seuil correspond au 
premier cycle de la phase d'efficacite decroissante pour 
lequel I'efficacite est inferieure au seuil d'efficacite 
minimal tol§r§ (par exemple O^l.Eo). 



Plus g6n6ralement, la valeur du seuil d'efficacit§ peut 
pr§ferentiellement etre comprise entre 0 et 0,5, et une 
efficacit6 de PGR inferieure a cette valeur est 
potentiellement biais6e par des ph6nom6nes d' inhibition 
incontrdl6s . 

Dans !• exemple repr6sent§ sur la figure 7, la valeur seuil 
d'efficacit6 a §te fixee a 0,02 (soit 2% de Eq) . Le cycle 
seuil Ca correspondrait au cycle n=36. On a repr6sente sur 
cette figure 7 les mesures de fluorescence 6mises 
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ajustees. On constate une correlation satisf aisante du 
module (trait plein) avec les mesures expferimentales 
(marquees "o") jusqu'au cycle n=36- La figure 8 montre 
aussi una bonne correlation avec les mesures 
exp6rimentales des predictions de 1 ' ef f icacite, obtenue A 
partir de la figure 7 et de la mod61isation de la 
fluorescence mesuree et ajustee." . • • - 

Les 6tapes principales de ce mode de realisation peuvent 
se resumer comme suit, en reference S la figure 9. 

A I'etape de debut 70, on a obtenu les valeurs de 
quantites de fluorescence mesurees et ajustees par rapport 
au bruit de fond en fonction du nombre de cycles n, comme 
represente sur la figure 6A. 

A I'etape 71, on calcule une approximation de I'efficacite 
de la reaction au cycle n grace k la formule (2) ci-avant, 
pour tout cycle n = 1, 2, (N-1) . 

A I'etape 72, on determine le cycle Cg minimal de fagon a 
ce que les deux conditions suivantes sont satisfaites : 

- le cycle Cs se situe dans la phase d'efficacite 
decroissante, 

- I'efficacite au cycle seuil est inferieure k la valeur 
seuil d'efficacite Es (par exemple Eg = 0,1. Eo) : 

Ecs ^ Es 

On peut deja eliminer les points NEG pour lesquels . 
I'efficacite est inferieure ^ Es. 
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A l'6tape 73, on modSlise la courbe de fluorescence 6mise 
ajust6e dans la d6croissance de 1 'ef f icacit6, sur la 
portion de cycles Ceep = (Cg - 5) ^ Cg, grSce ^ I'expression 
(8) dans laquelle on a fait I'hypothdse sur la 



Ensuite, le test 74 sur la valeur estim6e E'^ de E'o et la 
decrementation a l'6tape 75 de la valeur du cycle de 
basculement Ceep visent S diminuer la valeur recherch^e de 
Cbep d'un pas P (qui peut §tre §gal A . 1) et A r6p6ter 
l'6tape 73 tant que la valeur de E'^ est inf6rieure S 1. 

Ensuite, lorsque la valeur d'efficacite estim6e d6passe la 
valeur 1 (fl§che n en sortie du test 74), on incremente la 
valeur de I'indice Ceep d'un pas h (qui peut §tre une 
fraction d'unit6) & I'etape 76 et on proc6de ^ une 6tape 
77 de modelisation de la fluorescence Fn du type de celle 
men^e A l'6tape 73, Tant que I'efficacite estim^e E'o est 
sup6rieure ou 6gale a 1 au test 78, on rep6te les etapes 
7 6 A 78. Lorsque I'efficacite estimee prend une valeur 
inf6rieure 1, on conserve les paramdtres estim6s 

(Fo,Fo.P'>Ceep) a I'etape de fin 79. 

A cette etape, on a obtenu finalement une valeur F,, 
seulement representative de I'effectif initial No dans 
I'echantillon d'int6r§t. On peut alors avoir recours ^ au 
moins un echantillon standard d'effectif initial Nq®^ 



5 



compensation e telle que 8=F'o : 
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connu, pour determiner, a I'etape 80^ I'effectif initial 
No dans I'^chantillon d'int6r§t. 

A cet effet, on obtient une valeur de mesure d'un effectif 
5 initial Fost dans un echantillon standard d' effectif 
initial connu Nost- On determine ensuite la valeur de 
I'effectif initial dans 1 ' 6chantillon d'int6x§t No par 
application d'une m§me relation de proportionnalit6 entre 
I'effectif initial No et sa mesure F'of dans 1 ' 6chantillon 
10 standard d'une part, et dans 1 ' Echantillon d'int6r§t 
d' autre part. 

En d'autres termes, a I'Stape 80 de la figure 9, on peut 
determiner la valeur No de I'effectif initial dans 
15 1' echantillon d'int6r§t par une simple regie de trois du 

type : No^F'o (n 

ost'^'ost)' impliquant I'effectif initial dans 

le standard Nost et un rapport des fluorescences corrig6es, 
compens6es et estim6es par ajustement du module de 
fluorescence, pour 1' echantillon d' interet et pour 
20 1' echantillon standard. 

On comprendra ainsi qu'un seul standard devrait etre 
suffisant pour determiner I'effectif initial de la 
population d' interet dans 1 ' echantillon d' interet, selon 
25 un avantage que procure 1' invention. 

Neanmoins,' dans une variarite qui peut aussi etre pr6vue si 
necessaire, on obtient les valeurs respectives de mesures 

d'effectifs initiaux F'qj, dans plusieurs echantillons 
30 standard d'effectifs initiaux connus Noat. On etablit 
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ensuite une d6piendance entre les effectifs initiaux Noat 
des echantillons standard et leurs valeurs respectives de 
mesure d'effectif initial F^^ • Ensuite, aprds la 
d6termination de la valeur de mesure de I'effectif initial 
5 pour l'§chantillon d'intergt , on determine I'effectif 
initial de la population d'inter§t No. par interpolation 
sur cette d6pendance . On comprendra que cette dependance 
peut etre encore typiquement une regression du type 
represents sur la figure 5, mais sur laquelle on reportera 
10 en ordonnees les fluorescences initiales Fq^, des standards 
et de 1' echantillon d' int§rSt F^ (ou leurs valeurs 
logarithmes respectives), plutot que leur valeur d' indice 
de basculement Cebp- 

15 D6s lors que I'on fait appel k un ou plusieurs standards, 
on peut prSvoir un ou plusieurs 6chantillons standard 
d' effectifs initiaux respectifs connus Nostf auxquels on 
applique la succession de reactions d'amplif ications, 
sensiblement dans les memes conditions que pour 

20 1 '6chantillon d'interet. Ensuite, la determination de leur 
valeur de mesure d'effectif initial F^^^ se fait par 
Qomparaison des variations param6tr§es avec les mesures 
expSrimentales, comme pour 1' echantillon d'int6r§t. 

25 En d'autres termes, on applique bien entendu les m§mes 
calculs pour les quantit6s de fluorescence mesurSes et 
ajustSes sur le ou les standard (s) et pour 1' echantillon 
d'int6r§t. On estime la quantite de fluorescence F^, avant 
le premier cycle, pour le ou les standard (s) selon le m§me 
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procede de determination de Fq pour I'echantillon 
d' interet , decrit ci-avant . 

Un troisieme mode de realisation, correspondant a une 
variante du second mode de realisation ci-avant consists 
globalement a ajuster le modele de I'efficacite En par 
rapport aux mesures exp^rimentales et ^ injecter ensuite 
ce module d' ef f icacite a juste dans le modele de I'effectif 
courant Nn ou encore dans le module des mesures Fn. Ce 
troisieme mode de rfealisation se r§sume corame suit. 

La variation parametr6e construite ^ I'etape b) est 
representative du rendement, et, ^ I'fetape c) , on 
determine une variation exp6rimentale du rendement a 
partir des mesures exp.6rimentales, pour comparer la 
variation parametr6e A la variation expferimentale. 
Ensuite, pour acceder a un param6tre repr6sentatif de 
I'effectif initial No, ^ I'etape d) : 

dl) on determine une seconde variation param6tr6e, 
representative de I'effectif courant Nn dans 
I'echantillon d' interet, et mettant en jeu au moins le 
param^tre repr^sentant I'indice de basculement Cbep/ 
ainsi qu'un parametre repr6sentatif de I'effectif 
initial No, 

d2) on applique k cette seconde variation une valeur du 
parametre de I'indice de basculement Ceep d6termin6 k 
1 'etape c) , 

d3) et on a juste au moins le paramdtre representatif de 
I'effectif initial No par comparai'son directe de la 
seconde variation avec les mesures exp6rimentales - 
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On peut avantageusement : 

- A l'6tape d2), appliquer une valeur grossiere de 
I'indice de basculement Ceep. comme on a decrit sa 
detection en reference k la figure 4B ci-avant, tandis 
que 

- ^ l'6tape d3), on peut af finer ensuite la valeur de 
I'indice conjointement avec I'ajuatement du pBramStre 
representatif de I'effectif initial No. 

Finalement, on retiendra simplement que le second mode de 
realisation des figures 7 et 8, actuellement pr6f§r6, se 
distingue de ce troisidme mode de realisation par le fait 
que I'on ne cherche pas A mener les correlations sur 
I'efficacite, mais on utilise simplement le module 
mathematique de variation de 1 • ef f icacite pour modeliser 
et af finer 1' estimation de la fluorescence corrig6e et 
compensee . 

Bien entendu, la pr^sente invention ne se limite pas a la 
forme de realisation decrite ci-avant ^ titre d'exemple ; 
elle s'etend k d'autres variantes . 

Ainsi, on comprendra que la presente invention peut 
s' appliquer en outre a la quantification relative, 
notamment par PCR. Dans cette application, on prevoit 
d' amplifier, outre la population d'interet, une population 
de reference soit conjointement dans le meme support, soit 
separement. On relive respectivement : 

- les mesures expirimentales representatives de 
I'effectif de la population d'int6r@t, et 
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- les mesures exp6rimentales representatives de 
I'effectif de la population de r6f6rence. 

Le proc6de peut se poursuivre alors par 1 • application^ a 
5 la population de reference/ des §tapes a), b) et c) , 
tandis que I'etape d) consiste siinplement ^ determiner un 
ratio- entre les* effectifs" iTiitiaux' respectifs -de ' la 
population d'int6ret et de la population de r6f6rence, 

10 On indique que la quantification * relative peut etre 
utilis6e pour analyser 1' expression d'un gSne d'interet au 
cours du dSveloppement d'un organisme. Pour corriger 
notamment les variations de quantity et de qualite entre 
les echantillons preleves sur 1' organisme ^ differents 
instants, on analyse, en plus du gene cible d'interet, un 
gdne de reference connu pour avoir un niveau d' expression 
stable au cours du d§veloppement • 

II s'agira finalement de comparer les ratios entre 
les differents 6chantillons preleves. 



15 



20 



25 



Pour aboutir aux resultats souhait6s, deux strategies sont 
possibles. 

Celle de I'art anterieur basee sur la detection du cycle 
Ct se d6roule habituellement comme suit. On determine, 
pour chaque echantillon pr^leve ^ differents instants tO, 

N 

tl, t2, tn, le ratio -rp^ # en utilisant au moins un 

standard (c' est-^-dire un Echantillon dont No^tte et N^^^ 
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sont connus), ce qui revient a faire deux quantifications 
absolues successives suivies d'un calcul de ratio. 

Une autre strategic particuli6rement avantageuse dans, le 
contexte de 1' invention consiste k determiner, pour chaque 
6chantiHon pr61ev6 diff6rents instants tO, tl, t2,..., tn 

directement le ratio .[^"^"^ grace ^ la formule : 



(n 1 

V ^"'^ J id. ,0 



Dans cette seconde realisation qui n' utilise finalement 
que le param^tre Fo, en combinant la technique de 
1' invention, aucun ^chantlllon standard n'est n§cessaire, 
de fagon particulidrement avantageuse. 

On se refere maintenant ^ la figure 10 sur laquelle on a 
represents une installation pour la mise en ceuvre du 
proc6d6 selon 1' invention. Elle comporte un support SUPP 
comportant ici un puits comprenant 1 'echantillon d'interet 
ECH et un puits comprenant par exemple un Schantillon 
standard not6 St. Le support SUPP est enferme dans une 
enceinte ENC par exemple iquipee de moyens de chauffage 
(non repr6sent6s) pour appliquer au standard et A 
l'6chantillon d'int6ret une reaction de PGR. 
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Dans I'exemple decrit, on pr6voit pr6f 6rentiellement 
d'appliquer des mesures de quantit6s de fluorescence 
§inises ^ chaque cycle par le standard St et par 
l'6chantillon d' inters t ECH, A cet effet, on introduit 
dans les puits un r^actif choisi et on 6claire les 
echantillons ^ partir d'une lampe (par exemple halogene- 
tungstSne) • dans le but - de ■ mesurer des -quantites 
respectives de fluorescence 6inanant de 1 '^chantillon 
d'interet et de 1' 6chantillon standard, ^ chaque cycle de 
PGR qui sera appliqu§. En outre, un appareil de d6tection 
de fluorescence comporte par .exemple un objectif 11 pour 
collecter la lumiere issue de la fluorescence, ainsi que 
des moyens de comptage de photons 10, par exemple une 
camera CCD, et/ou des photomultiplicateurs, pour compter 
15 la fluorescence emise ^ chaque cycle de PGR, par 
I'echantillon d'interet et le standard. Ainsi, la 
fluorescence emise par chaque puits est avantageusement 
f ocalisee par 1 'objectif 11, puis d6tectee 
preferentiellement par une camera CCD 10 connectee ^ une 
carte d ' acquisition 21, par exemple de type PCMCIA prevue 
dans une unite centrale 20 d'un ordinateur. 

Cet ordinateur est alors relie aux moyens de comptage 10 
pr6cites pour recevoir des signaux representatif s de 
mesures de quantite de fluorescence detectee a chaque 
cycle de PGR et traiter ces signaux pour determiner un 
effectif initial de la population d'interet avant le 
premier cycle, par la mise en oeuvre du precede au sens de 
1 • invention. 



20 



25 



30 
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On indique que l'unit6 de traitement comport e 
typiquement : 

- une carte d' acquisition 21 connect6e aux moyens de 
comptage 10/ 

- une mSmoire de travail 25 (par exemple de type RAM) 
pour le stockage temporaire et le traitement des 
signaux pr6cit6s, ... „ 

- une m6moire permanente 24 pour le stockage du produit 
programme d'ordinateur au sens de 1 ' invention et des 
donnees traitees et pretes ^ exploitation, par exemple 
pour un diagnostic ulterieur, 

- le cas 6cheant, un lecteur 22 d'un support memoire tel 
qu'un CD-ROM, ou autre, sur lequel peut Stre stocke 
initialement le produit programme d'ordinateur, 

15 - eventuellement, une interface de communication 26 vers 
un ssite local ou distant (liaison 28), par exemple pour 
la transmission des donnees traitees pour un diagnostic 
sur le patient realise a distance, 

- une interface graphique 27 reliee typiquement ^ un 
20 moniteur d'affichage 30, 

- et un processeur 23 pour gerer les interactions entre 
ces differents ^quipements. 

On indique que I'ordinateur peut comporter aussi des 
25 organes de saisie tels qu'un clavier 41 et/ou une souris 
42 relies k 1' unite centrale 20. • 

On retlendra toutefois que 1' installation au sens de 
1' invention coraporte globalement : 
30 - un support SUPP d» echantillon au moins pour 
1' 6chantillon d'int6r§t. 
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- un premier appareil ENC pour appliquer, au moins d la 
population ci'inter§t dans 1 ' 6chantillon d'intSrSt, 
ladite succession de reactions d' amplification, 

- un second appareil 10 pour relever des mesures 
5 representatives de I'effectif courant de la population 

d'intdrSt, et 

-"de-s moyens inf ormatiques" 20 propres A r6cup6rer des 
signaux de mesures du second appareil 10 et S mettre en 
oeuvre tout ou partie des etapes du precede au sens de 
10 1' invention. 

A cet effet, un produit programme d'ordinateur peut animer 
ces moyens inf ormatiques . Ce programme peut etre stocke 
dans und memoire de 1' unite de traitement 20, ou sur un 

15 support memoire amovible (CD-ROM ou autre) et destine k 
coop6rer avec le lecteur de I'unitfe de traitement. Le 
programme d'ordinateur au sens de 1' invention comporte 
alors des instructions pour la mise en oeuvre de tout ou 
partie des §tapes du proc^d^ selon 1' invention. Par 

20 exemple, I'algorithme de ce programme peut etre represents 
par un organigramme Equivalent au schema de la figure 9. 
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Revendications 



1. Proc6d6 mis en oeuvre par des moyens informatiques pour 
quantifier, de fagon absolue et/ou relative, une 
population initials d'acides nucl6iques dans un 
6chantillon d'int6ret soumis ^ une succession 
d^applications • d'une -reaction • d- amplification- de la 
population, au cours desquelles on relive des mesures 
exp§rimentales representatives d'un effectif courant de la 
population au moins dans 1 • 6chantillon d'int6ret, 

caract6rise en ce qu'il comporte les 6tapes : 
a) prevoir un modele de rendement (E„) de la reaction 
d- amplification en fonction de la succession 
d' amplifications, ce module comportant : 

° Ph^se sensiblement constante (Eo) pour une 

premiere partie des applications de la reaction, 
et 

o une phase non constante pour une seconde partie 
des applications de la reaction, 
les premiere et seconde parties etant separ6es par une 
region de basculement du rendement entre les phases 
constante et non constante, incluant un indice 
d- application (Cbep) correspondant sensiblement au 
basculement, 

b) exprimer, en utilisant le module de rendement, une 
relation faisant intervenir au moins 1- indice de 
basculement et un paramStre representant 1' effectif 
initial dans 1 • echantillon d'int6rSt (No;Fo), 

c) par comparaison avec les mesures exp6rimentales, 
determiner au moins 1' indice de basculement, et, dans 



30 
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une 6tape d) ulterieure ou conjointe, en d^duire 
I'effectif initial dans 1 '^chantillon d'int6rSt (No). 

2. Proc6d6 selon la revendication 1, caract6ris6 en ce 
5 que : 

a I'etape b) , on exprime, en utilisant le module de 
rendement, une variation -param§tr6e en -fonction de 
ladite succession d' amplifications, mettant en jeu au 
moins un param^tre repr6sentant I'indice de basculement 
10 (Cbep) , et 

A I'etape c) , par comparaison avec lesdites mesures 
experimentales, on determine au moins ce param^tre 
repr^sentant I'indice de basculement (Cebp) . 

15 3. Procede selon la revendication 2, caracteris6 en ce que 
ladite variation param§tree est representative de 
I'effectif courant (N„) dans 1 ' 6chantillon d'int§ret, 
en ce que ladite variation met en jeu en outre un 
paramdtre representatif de I'effectif initial dans 

20 l'6chantillon d'interet (No), 

et en ce que, aux etapes c) et d) , on determine 
sensiblement conjointement les paramStres representatif s 
dudit indice d' application (Cebp) et de I'effectif initial 
(No) dans 1 ' echantillon d'int^rSt. 

25 

4. Procede selon la revendication 2, caract6ris§ en ce que 
ladite variation param6tr6e est representative du 
rendement, et en ce que,, a I'etape c) , on determine une 
variation exp6rimentale du rendement ^ partir desdites 
30 mesures experimentales, pour comparer la variation 
paramitree d la variation exp6rimentale . 
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5. Proc6de selon la revendication 4, caract§ris§ en ce 
que, i I'Stape d) : 

dl) on determine une seconde variation param6tr§e, 
representative de I'effectif courant (Nn) dans 
I'echantillon d'int^rat, et mettant en jeu au molns le 
■ • paranmStre -reprfesentant - ledit indice d'-cipplrcatlon 
(Ceep) f alnsi qu'un param^tre repr6sentatif de 
I'effectif initial (No) dans I'echantillon d'int6r€t, 

d2) on applique ^ la seconde variation une valeur du 
paramdtre de 1' indice (Ceep) determine S l*§tape c) , 

d3) et on a juste au moins le parametre representatif de 
I'effectif initial (No) par comparaison directe de la 
seconde variation avec les mesures experimentales • 

6. Precede selon la revendication s/ caract^rise en ce 
que : 

- d I'etape d2) , on applique une valeur grossiere de 
1' indice (Ceep) ^ tandis que 

- a I'etape d3) , on affine la valeur de 1 ' indice 
conjointement avec 1 ' a justement du parametre 
representatif de I'effectif initial (No) • 

7. Precede selon la revendication 3 ou la revendication 6, 
caracterise en ce que ladite variation parametr^e ou, 
respectivement, ladite seconde variation parametree : 

- est representative desdites mesures exp6rimentales, et 

- comporte un paramdtre correspondant ^ une valeur de 
mesure (Fq) representative de I'effectif initial/ 

et en ce que I'on determine la valeur de mesure de 
I'effectif initial (Fo) par comparaison desdites 
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variations . parametr6es (Fn) avec les mesures 
exp^rimentales . 

8. Proc6d6 selon I'une des revendications pr^c^dentes, 
5 caract6ris6 en ce qu'on applique une 6tape prealable de 

traitement des mesures exp§rimentales comprenant une 
s-ous-traction d'un • bruit- -de fond de - mesures - et 
d' introduction d'une compensation pour tenir compte d'une 
mesure non nulle (e) representative de I'effectif initial. 

10 

9. Precede selon I'une des revendications 7 et 8, 
caracterise en ce qu'on obtient une valeur de mesure d'un 
effectif initial (Fq^^) dans un 6chantillon standard 
d'effectif initial connu (No^^) , 

15 et en ce qu'on determine la valeur de I'effectif initial 
dans 1 'echantillon d'interet (No) par application d'une 
meme relation de proportionnalite entre I'effectif initial 
et sa mesure, dans 1 * echantillon standard d'une part, et 
dans 1' echantillon d'interet d' autre part. 

20 

10. Procede selon I'une des revendications 7 et 8, 
caracterise en ce qu'on obtient des valeurs respectives de 
mesures d'effectifs initiaux (Fq^^) dans des echantillons 
standards d'effectifs initiaux connus (Nq^*") , 

25 et en ce que : 

- on etablit une dfependance entre les effectifs initiaux 
(No^^) des 6chantillons standard et leurs valeurs 
respectives de mesure d' effectif initial (Fo^^) r 

- et, apres la determination de la valeur de mesure de 
30 I'effectif initial pour 1 ' Echantillon d'inter§t (Fo) , on 
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d6termine I'effectif initial de la population d'inter§t 
(No) par interpolation sur ladite dependance. 

11. Proc6d6 selon I'une des revendications 9 et 10, 
caracteris6 en ce que I'on prSvoit un ou une pluralite 
d ' 6chantillons standard d ' ef f ectif s initiaux respectif s 
connus - (No'**') , • auxquel's -on- applique ladite - succession 
d' amplifications, sensiblement dans les mimes conditions 
que pour 1 • echantillon d'interet, ainsi que la 
determination de leur valeur de mesure d'effectif initial 
(Fo^^) par comparaison des variations param6trees (Fn''*') 
avec les mesures exp6rimentales . 

12. Procede selon la revendication 3, caract6rise en ce 
que I'on prevoit une pluralite d • echantillons standard 
d'effectifs initiaux respectif s connus, auxquels on 
applique ladite succession d' amplifications, sensiblement 
dans les memes conditions que pour 1 ' echantillon 
d'intferet, ainsi que la determination de leurs indices 
respectif s selon les etapes a), b) et c) , 

et en ce que, a I'^tape d) : 

- on etablit une dependance entre les effectifs initiaux 
des echantillons standard (No^^) et leurs indices 
(Ceep^^) / 

25 - et, aprds la determination de I'indice pour 
• l'6chantillon d'interSt (Cebp) , on determine I'effectif 
initial de la population d'int^ret (No) par 
interpolation sur ladite dependance • 



15 



20 



30 



13. Proc6d6 selon I'une des revendications 1^8, 
caract6ris6 en ce que, pour une quantification relative. 
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on pr6voit, outre la population d'int6ret, une population 
de r§f§rence soumise ^ une succession d' applications de 
reaction d' amplification, et on relive respectivement : 

- les mesures experimentales representatives de 
5 I'effectif de la population d'interet, et 

- les mesures exp6rimentales representatives de 
" I'effectif de la populatrion de reference, 

le procSdfe se poursuivant par 1 ' application, d la 
population de r6f6rence, des etapes a), b) et c) , tandis 
10 que l'6tape d) consiste ^ determiner un ratio entre les 
effectifs initiaux respectifs de la population d'int^ret 
et de la population de refe'rence. 

14, Proc6d6 selon I'une des revendications pr^cedentes, 
15 caracterise en ce que : 

- on exprime les mesures experimentales sous la forme 
d'une variation experimentale du rendement en fonction 
de ladite succession d ' applications, 

- on obtient une variation experimentale du rendement en 
20 fonction de ladite succession d ' applications, 

comportant : 

* une premiere region sensiblement bruitee pour des 
indices d * applications faibles, 

* suivie d'une seconde region moins bruitee pour des 
25 indices d * applications plus eieves, 

15 • Precede selon la revendication 14, dans lequel ladite 
phase non constante du rendement est decroissante, 
caracterise en ce que : 
30 - on estime une valeur grossiere de phase constante du 
rendement (Eo) , et 
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- au moins pour la recherche de I'indice (Cbep) dans 
ladite region de basculement, on ignore une partie au 
moins des mesures dans ladite seconde region moins 
bruitee et dont le rendement estime est inferieur a une 
5 valeur seuil, de pr6f6rence k une fraction de la phase 

constante (Eo) . 

16. Proc6d6 selon I'une des revendications 14 et- 15, dans 
lequel ladite phase non constante du rendement est 

10 decroissante, caract6ris6 en ce que I'on repdre ladite 
region de basculement, dans le sens des indices 
d' applications decroissants ^ partir de ladite seconde 
region moins bruitee, en d6tectant un indice grossier 
d' application d" amplification dont le rendement depasse 

15 sensiblement une valeur predetermin6e (Eo=l) . 

17. Proc6d6 selon la revendication 16, caract6rise en ce 
qu'on affine 1' estimation de la valeur, le cas 6cheant 
f ractionnaire, dudit indice d' application (Ceep) dans la 

20 region de basculement, dans le sens des indices 
d ' applications croissants ^ partir de I'indice grossier, 
en detectant un indice d * application d ' amplification dont 
le rendement vaut environ , ladite valeur predeterminee 
(Eo«l). 

25 

18. Proced^ selon I'une des revendications precedentes, 
•caracterise en ce que ladite phase non constante du 

rendement est modelisee par une exponentielle d^croissante 
comportant un param^tre de decroissance (P), et en ce que 
30 ledit parametre de decroissance (P) est determine, a 
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I'etape c) avec I'indice dans la region de basculement 
(Ceep) f par comparaison avec las mesures exp6rimentales . 

19. Proc6d6 salon I'une des revendications pr6cedentes, 
5 dans lequel la reaction d' amplification est una reaction 
de PGR {"Polymerase Chain Reaction") en temps r6el. 



20. Proc6d6 salon I'une des revendications pr6c6dentes, 
dans lequal lesdites mesures sont des quantitfes de 

10 fluorescence 6mises, 

21, Installation pour la mise en cBuvre du precede selon 
I'une des revendications pr6c6dentesr caract6ris6e en ce 
qu'elle comporte : 

15 - un support d' echantillon au moins pour 1 ' echantillon 
d'inter§t, 

- un premier appareil pour appliquer, au moins ^ la 
population d'int^ret dans 1 • Echantillon d'interet, 
ladite succession de reactions d' amplification, 

20 - un second appareil pour relever des mesures 
representatives de I'effectif courant de la population 
d'interet, et 

- des moyens inf ormatiques propres a recup6rer des 
signaux de mesures du second appareil et a mettre en 

25 ceuvre tout ou partie des etapes du precede selon I'une 

des revendications precedentes. 

22. Produit programme d' ordinateur , destin6 a etre stock6 
dans une memoire d' une unite de traitement, ou sur un 
30 support memoire amovible at destine ^ coop6rer avec un 
lecteur de ladite unite de traitement, caract6ris§ en ce 
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qu'il comporte des instructions pour la mise en oeuvre de 
tout ou partie des etapes du proc6d§ selon I'une des 
revendications 1 a 20. 



! 
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ABREG6 



Procede, ina-tallat:xon et programiae d'ordinateur pour 
es timer I'effecbif initial d'mie population d'acidea 
nuclei<me3^ notamment par PGR 

Pour estimer une population initials d'acides" nucleiques 
dans un 6chantillon d'intferSt, notanunent par PGR : 

a) on pr6voit un module d'efficacite (£„) de PGR 
comportant une phase constants (Eo) suivie d'une phase 
non constants, s6par6es par une region de basculement 
incluant un indice de basculement (Cbep) r 

b) on exprime, en utilisant le modele d'efficacite, une 
relation entre 1' indice de basculement et un parametre 
representant I'effectif initial, 

c) par comparaison avec les mesures exp6rimentales, on 
determine 1' indice de basculement, et on en d6duit 
I'effectif initial dans 1 ' echantillon d'interet, 

(FIGURE 4B) 



Coov orovided by USPTO from the IFW Imajie Database on 01/19/2006 



■is 





4 ,st4 /St^- rat St3 

(a .a) 




CH Cirz Ctinir Cf:3 




BEST AVAILABLE COPV 



('u .£1,) 



CoDv provided by USPTO from the IFW Imaae Database on 01/19/2006 




CoDV provided by USPTO from Ih© IFW Image Database on 01/19/2006 




BEST AVAILABLE COP' 



Copy provided by USPTO from the IFW Image Database on 01/19/2006 



Us 




pr\'G -3 



75 



1 



Mod(F„) 




Dec(CeEp); P 


^ 




78 












^ o 






n 




savo (F'q, 1 










^80 




No 





,79 



BiEST AVAIUBLE COPY 



Copy provided by USPTO from the IFW Image Database on 01/19/2006 




^lEST AVAU-ABLE COP>^ 



CoDv Drovided bv USPTO from the IFW Imaqe Database on 01/19/2006 



Document made available under the 
Patent Cooperation Treaty (PCT) 

International application number: PCT/EP2005/0 13021 
International filing date: 18 November 2005 (18.11.2005) 

Document type: Certified copy of priority document 

Document details: Country/Office: FR 

Number: 0412471 

Filing date; 24 November 2004 (24.1 1.2004) 

Date of receipt at the International Bureau: 07 March 2006 (07.03.2006) 

Remark: Priority document submitted or transmitted to the International Bureau in 
compliance with Rule 17.1(a) or (b) 




World Intellectual Property Organization (WIPO) - Geneva, Switzerland 
Organisation Mondiale de la Propriete Intellectuelle (OMPI) - Geneve, Suisse 



